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Approximations made to estimate two-nucleon transfer probabilities in ground-state to ground- 
state transitions and physical interpretation of these probabilities are discussed. Probabilities are 
often calculated by approximating both ground states, of the initial nucleus A and of the final nucleus 
yl±2 by the same quasiparticle vacuum. We analyze two improvements of this approach. First, the 
effect of using two different ground states with average numbers of particles A and A±2 is quantified. 
Second, by using projection techniques, the role of particle number restoration is analyzed. Our 
analysis shows that the improved treatment plays a role close to magicity, leading to an enhancement 
of the pair-transfer probability. In mid-shell regions, part of the error made by approximating the 
initial and final ground states by a single vacuum is compensated by projecting onto good particle 
number. Surface effects are analyzed by using pairing interactions with a different volume-to-surface 
mixing. Finally, a simple expression of the pair-transfer probability is given in terms of occupation 
probabilities in the canonical basis. We show that, in the canonical basis formulation, surface 
effects which are visible in the transfer probability are related to the fragmentation of single-particle 
occupancies close to the Fermi energy. This provides a complementary interpretation with respect 
to the standard quasiparticle representation where surface effects are generated by the integrated 
radial profiles of the contributing wave functions. 

PACS numbers: 21.10.Pc, 21.10.Ro, 21.60.Jz, 25.40.Hs, 27.60.+j 



I. INTRODUCTION 

The link between Cooper-pair superfluidity in nuclei 
and cross sections associated to pair-transfer reactions 
has been extensively discussed in the literature starting 
from the early work of Broglia and collaborators in the 
70s In particular, the relation between the character- 
istics of the pairing correlations and the transfer prob- 
abilities in two-particle transfer reactions has been an- 
alyzed 0. Recently, there is a renewal of interest on 
experimental 1-nucleon (In), 2-nucleon (2n) and more 
generally multi-nucleon transfer channels at bombarding 
energies close to the Coulomb barrier [s*]. From the the- 
oretical side, in the last decade, new microscopic calcula- 
tions have been developed, essentially in the framework 
of the Hartree-Fock-Bogoliubov (HFB) -f quasiparticle 
random-phase approximation (QRPA) theory to inves- 
tigate the properties of the 0+ [4] and 2+ Q excita- 
tion modes associated to 2n addition or removal during 
transfer reactions. A similar analysis has been performed 
also within the time-dependent HFB model in the small- 
amplitude limit 0. The possibility to use pairing vi- 
brations to constrain the pairing interaction employed 
in HFB-based calculations has been explored recently 
0, The objective of these last studies was to sug- 
gest an experimental measurement adapted to identify 
the surf ace/ volume mixing character of the pairing in- 
teraction. Two-neutron (p, t) transfer reactions in very 
neutron-rich Sn isotopes have been indicated as a good 
experimental candidate. 

In the recent Ref. @, the effects of the surface/volume 



nature of the pairing interaction on 2n transfer have 
been extensively analyzed. The enhancement of transfer 
probabilities at the surface has been predicted (when a 
surface-peaked interaction is used) in the transition from 
the ground state (GS) of the nucleus with mass A to the 
GS of the nuclei with masses A ± 2 in Sn isotopes be- 
yond = 82. Similarly to what done in other recent 
estimates of the pair-transfer probability using the mi- 
croscopic HFB approach Q, the strength associated to 
these transitions has been calculated in Ref. Q with an 
approximate formula where only the wave functions of 
the nucleus A enter. This approximate treatment differs 
from the original formulation given in Ref. jjy] where the 
components of both nuclei A and Azt2 appear. Accord- 
ingly, in other recent works not based on HFB ^To'], in 
the expressions of the two-particle transfer spectroscopic 
amplitudes the wave functions of the two nuclei A and 
A + 2 appear It is worth mentioning that, even the 
formulas given in Ref. [l[ are approximate expressions 
and, as far as we know, the underlying approximation 
needs to be clarified. 

One goal of the present work is to discuss different level 
of approximation used to estimate two-nucleon transition 
probabilities. In this work we derive the expressions to 
be used for the two-particle GS GS transfer probabil- 
ities in the framework of the HFB model. We compare 
our results with those obtained with the model of Shi- 
moyama and Matsuo @ and analyze the differences in 
mid-shells and at shell closures. Following their work, 
different values for the surface/ volume mixing parame- 
ter are used in the pairing interaction and application is 
made for the chain of Sn isotopes. By using the canoni- 
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cal basis representation of the HFB model we provide a 
complementary interpretation of surface effects. Finally, 
we quantify the effect of particle number restoration on 
two-neutron transfer probabilities. 

The article is organized as follows. In Sec. II the gen- 
eral scheme of the present calculations is presented. In 
Sec. Ill, an expression for the transfer probability is de- 
rived in an approximated framework analogous to that of 
Ref. in the quasiparticle (coordinate) (Subsec. Ill- A) 
and in the canonical basis (Subsec. III-B) representation. 
The interpretation of the pair-transfer probability is dis- 
cussed in Subsec. III-C. In Subsec. III-D the correspond- 
ing results are shown and compared. Improved formulas 
for the probability are derived in both the quasiparticle 
(Subsec. IV-A) and the canonical basis (Subsec. IV-B) 
formulations. Results are presented and commented in 
Subsec. IV-C. Finally, projection techniques are applied 
in the canonical basis case. The effect of particle num- 
ber conservation on two-neutron transition probabilities 
is discussed (Sec. V). In Sec. VI, a summary is drawn 
and perspectives are outlined. 



fm ^ for a mixed interaction {rj = 0.65). The Skyrme in- 
teraction which is employed in the present calculations 
is SLy4 
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The HFB calculations are performed in 
coordinate representation with a box discretization and 
a box radius equal to 20 fm. The pair-transfer proba- 
bilities are derived both in the quasiparticle framework 
and in the canonical basis. In the latter case, it will be 
shown that the integrals of the radial parts of the wave 
functions disappear in the expressions of the pair-transfer 
probabilities (which depend only on the occupation num- 
bers). This allows us to check whether the differences 
found with different pairing potential profiles (enhance- 
ment effects when a surface-peaked interaction is used) 
are artificial effects which are provided only by the for- 
mulations where integrals of wave functions are present 
or genuine physical effects that are found in all types of 
formulations. 

Besides the approximation made on the many-body 
wave function within the HFB approach, further approx- 
imations have to be done to estimate the pair-transfer 
probability. These approximations are analyzed in the 
following sections. 



II. FRAMEWORK OF THE PRESENT 
CALCULATIONS 

The removal and addition GS — > GS pair-transfer am- 
plitudes are written as. 



-'gs 



(GSA_2|*,'(ri,-(7i)*,(r2,a2)|GSA), (1) 



and 



rriAdd 



(2) 



respectively, where a and q represent spin and isospin, 
respectively. These matrices are assumed diagonal in 
isospin, that is q' = q. The states |GSa) and |GSyi±2) 
correspond a priori to the highly correlated ground states 
in the entrance and exit channel of the transfer reaction. 
For medium and heavy nuclei, only approximate wave 
functions can be obtained 

Here, the framework of the Skyrme-HFB model with 
a zero-range density-dependent pairing interaction Vpair 
is used with 



Vpair{ri,r2) = Vo 



1 — 77 



p(R) 

Po 



S{r,-V2), (3) 



where R = (ri+r2)/2. The parameters a and po are cho- 
sen equal to 1 and 0.16 fm^^, respectively. With a cutoff 
equal to 60 MeV in quasiparticle energies and a maximum 
value of j equal to 15/2, the parameter Vq is adjusted to 
reproduce the two-neutron separation energies of Sn iso- 
topes as already done in Refs. {L,^. Two different pair- 
ing interactions are employed in terms of surface / volume 
mixing to check the sensitivity of the results on two dif- 
ferent pairing radial localizations. The values of the pa- 
rameter Vq which are used here are Vq = -670 MeV fm~^ 
for a pure surface interaction (77 = 1) and Vq = -390 MeV 



III. PAIR-TRANSFER PROBABILITIES 
WITHIN AN APPROXIMATED SCHEME IN 
THE HFB THEORY 

A. Quasiparticle formulation 

In Ref. the removal and addition amplitudes are 
calculated by assuming that the ground state of the nu- 
cleus A±2 can be approximated by the ground state of 
the nucleus A in the evaluation of Eqs. ^ and By 
taking ri = r2 = r and by assuming spherical symmetry, 
this approximation leads to 



-^GS 



nlj 

where u'^{r) and v'^{r) are the radial parts of the upper 
and lower components of the quasiparticle wave functions 
for the nucleus A, respectively, and k is the anomalous 
density of the nucleus A; the index n runs over the num- 
ber of states in each {j,l) channel. To derive the above 
expression in terms of the u and v components, the Bo- 
goliubov transformations have been used (after having 
done the approximation \GSa±2) ^ \GSa)), which can 
be written in the following way, 

n 

(5) 

One may expect that the approximation leading to Eq. 
(4) is reasonable for mid-shell nuclei. However, at shell 
closures, the ground states of the nuclei A and AzL2 are 
expected to be quite different and the anomalous density 
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is zero due to the collapse of static pairing correlations. 
It can thus be interesting to check especially in these 
regions the validity of such an approximation. 

The pair-transfer probability calculated by using the 
amplitude given by Eq. (|4]) is written as 



pRcm 
^GS 



{A) 



pAdd 
^GS 



{A) 



drY,{2j + l)ut^{r)vit,(.r) 

nlj 



(6) 



B. Canonical basis formulation 

An alternative formulation can be developed in the 
canonical basis formalism. Let us introduce the canon- 
ical states, denoted by {^Pi} associated to the set of 
creation/annihilation operators {a|,ai} where i contains 
all quantum numbers. The corresponding occupation 
number is denoted by n^, where A is the mass of the 
nucleus. Accounting for the convention for the time- 
reversed states used in Ref. \Vm we have 



dr{^,(r,t)*t(r,t)+</5.(r,i)*^(r,i)}, 
dv { (r, i)vl/t (r, t) - V* (r, t)*^ (r, ;) } , (7) 



where ^' are the operators introduced in the Bogoliubov 
transformations. These equation can eventually be in- 



verted to give 

vl/t(r,t) = 5]{a,V:(r,t) + al^|(r,t)}, (8) 

i>0 

vl/t(r,;) = ^{a,V:(r,i) + al^|(r4)}. (9) 

i>0 

In the canonical basis, the HFB quasiparticle ground 
state takes a BCS-like form, 



Gs^)=nK+^M4)io) 



(10) 



j>0 



where = ^\ — nf and = \fnf ■ In what follows, 
we consider the addition pair-transfer probability; the 
expression for the removal probability can be derived in 
the same way. We have: 



ij>0 

{aJ(p*(r,i)-l-a-V|(r,;)}|GSA). (11) 



Similarly to what is done in Ref. [9|] and in the previ- 
ous section, one can eventually assume that (GSa-i-2 | can 
be replaced by (GSa|. By taking advantage of the fact 
that the single-particle states are canonical and by using 
the properties relating time-reversed states, one finally 
deduces that 



Tht\^) = T«g'"(r)^-^(GS^|aIat|GS^)(|<^,(r,t)P + |^.(r,i)|^ 



i>0 



^^nf(l-nf)(|(p,(r,t)P + |^.(r,i)p) 



j>0 



In the special case considered here, where the nucleus is assumed spherical, we finally obtain 

I 



nlj 



(12) 



(13) 



where (j)nij stands for the radial part of the canonical ba- 
sis component. This expression provides an alternative 
form of the transition amplitude and its numerical esti- 
mate should exactly match with the one obtained with 
Eq. (4). Interestingly enough, in this specific basis, the 
radial dependence of the wave function completely dis- 
appears in the transition probability 



pAdd 



(A) 



pRcm 



(A) 



(14) 



or, equivalently, in spherical symmetry, 



E 

nlj 



(2j + l)\/< (!-<■) 



(15) 
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Comparison between quasiparticle and 
canonical formulation 



In Fig. [U the removal (or addition) probabilities ob- 
tained by using Eqs. (6) and (15) are compared for the 
Sn isotopic chain and the two different employed pair- 
ing interactions. For the surface case, the nucleus -'^'^^Sn 
corresponds to the drip line nucleus (the two-neutron sep- 
aration energy changes its sign going from ^'^^Sn to ^'^^Sn 
in this case). As can be seen in this figure, the two sets 



1400 
1200 
1000 



^ 800 



GO 600 
400 
200 




' 1 


■ 1 ■ 1 

Surface 


' 1 




Mixed 




:^ 




■ / - 

/ / 
/ / 
11 


- 'X 

If 

£_■ 1 


>. 

\ V 

V » 
\* ' 

V '> 

■ 1 1 JL 


1 w 

1 / 
/ 

1 



wonder whether this is an artificial result depending on 
the technical details of the HFB calculations, where in- 
tegrals of the components of the quasiparticle wave func- 
tions are done to evaluate the probability. However, in 
Eq. ([15]), the radial dependence of the canonical-state 
wave functions is integrated out and completely disap- 
pears. We can thus argue that the differences in the tran- 
sition probabilities are not artificial and certainly con- 
tain genuine physical effects. These effects that lead to a 
larger transfer probability for the case of a pure surface 
interaction with respect to the case of a mixed interac- 
tion could be seen (in the canonical basis formulation 
where the probability is expressed in terms of occupa- 
tion numbers) as due to a different neutron occupancy 
fragmentation around the Fermi energy. The neutron 
occupation numbers around the Fermi energy are more 
fragmented in the case of a surface pairing interaction; 
in this latter case the Fermi energy is also closer to zero 
meaning that the system is less bound. This indicates 
that the last occupied states are closer to the continuum. 
This is illustrated in the upper panel of Fig. [5] where the 
neutron Fermi energy A„ is displayed as a function of A 
for both interactions. To estimate in a systematic way 
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Figure 1: Comparison of the removal (addition) probability 
obtained by using Eq. ([6]) and for the mixed pairing case (blue 
solid line) and the pure surface case (red dashed line). The 
results obtained using Eq. (|15p are also shown for the mixed 
(filled circles) and pure surface (blue filled squares) case 

of results can almost be superposed one to the other, as 
it should be if a full basis was used in both formulations. 
The small difference is due to the fact that in both cases 
the actual calculation is made with a set of states be- 
low a certain cutoff. The cutoff used in the quasiparticle 
space case cannot be easily connected to the cutoff in 
canonical single-particle space. Therefore, the inevitable 
slightly different cutoff choices lead to the (very) small 
difference. 

The surface-peaked pairing interaction systematically 
provides a larger probability. We notice that the transfer 
probability is zero in this approximation for the magic 
nuclei ^""Sn and ^'^^Sn due to the absence of static pairing 
correlation in these cases. 
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D. Interpretation of pair-transfer probability 

The canonical basis formulation of the pair transfer 
gives an interesting new insight in the interpretation of 
the transfer probability. Indeed, the differences between 
the pair-transfer probabilities which are obtained with 
different pairing forces are usually interpreted as a con- 
sequence of the radial features of the quasiparticle wave 
functions that appear in the integral of Eq. (6). One may 



Figure 2: Upper panel: Evolution of the neutron Fermi energy 
for the tin isotopic chain for mixed (filled circles) and surface 
(filled squares) pairing case. Lower panel: evolution of the 
entropy for the same nuclei. 

the fragmentation of the single-particle occupancies, the 
single-particle entropy, defined as 

S = -A:B^(7i,lnn, + (l-n,:)ln(l-n,)), (16) 

i 
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is shown as a function of A in the lower panel of Fig. 
2. The more spreading of occupation number, the larger 
should be this quantity. A larger fragmentation of the 
occupation numbers implies a more diffuse Fermi surface 
and thus more important surface effects. 



basis representation. The two gap expressions are 
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Figure 3: Evolution of A for the tin isotopic chain for mixed 
(filled circles) and surface (filled squares) pairing case. 

This figure illustrates that the trend of the entropy is 
actually comparable to the trend of the strength. In both 
cases the results obtained with the surface-peaked pair- 
ing interaction are systematically larger than the values 
associated to a mixed interaction. 

It is worth spending some words about the trend of the 
square of the pairing gap as a function of A. In Ref. 
[9| the authors compare and the pair-transfer strength 
in Sn isotopes. They expect that these quantities should 
be proportional. By comparing their results, they actu- 
ally observe that the two quantities have the same trend 
but that they are not proportional: the surface enhance- 
ment effect beyond the isotope ^^°Sn (in their case the 
drip line for the surface interaction is shifted to heavier 
isotopes with respect to the present calculations) is more 
strongly visible in the pair-transfer probability than in 
A^. We believe that the pair-transfer strength and the 
square of the pairing gap cannot be directly compared. 
These two quantities cannot scale in the same way since 
they are obtained by integrating different functions. Fur- 
thermore, the comparison depends also on the adopted 
definition of the pairing gap which is not unique. In Ref. 
the expression 



A = 



J drK{r)A{r) 
J drK.{r) 



(17) 



is employed, where k is the anomalous density. Another 
definition is also currently used namely. 



A = 



/ drp{r)A{r) 
J drp{r) 



(18) 



where p is the particle density. 

To clarify the connection between A and A on one side 
and between these quantities and the pair-transfer prob- 
abilities on the other side, let us go back to the canonical 



A = 



~ ^ - n.)v/nj(l - n^QCj ^^^^ 

J2, ^/niil-Ui) 



where 



a, = / drVafMr^'lVjirW 



(20) 



We see that only if Cy = 1 for any couples of states 
around the Fermi energy, we have: 



A^ = A^ = P, 



Rom/ Add 



GS 



{A). 



(21) 



Such a condition is very unlikely. A and A mix in a 
different way radial effects and fragmentation of occupa- 
tion numbers showing again that these quantities cannot 
a priori directly be connected one to the other and to 
transfer properties. 

In the case of a mixed pairing interaction the two def- 
initions of A are expected to provide similar results be- 
cause the radial localizations of p{r) and K{r) are not 
very different. The case of a surface-peaked pairing in- 
teraction is however different. Due to the very different 
radial profiles of the anomalous and particle densities the 
two gaps may be quite different in this latter case. In 
particular, if the expression of Eq. (fT8)) is used, the pair- 
ing gap can have a very different behavior than the pair 
transition probability. This is shown in Fig. 3 where one 
sees that the surface-case A^ is lower than the volume- 
case A^ (inversed behavior with respect to the transfer 
probability) between the isotopes ^'^^Sn and ^-^^Sn. 



IV. IMPROVED TREATMENT OF PAIR 
TRANSFER 

To obtain the expression ^ or ([T5|) , it has been neces- 
sary to assume that \GSa) — \GSa±2) which is at vari- 
ance with the original prescription [1]. In this section, 
we discuss how to obtain, in the framework of the HFB 
model, a formulation similar to that advanced in Ref. [1| . 



A. Quasiparticle formulation 

We do not adopt the approximation used in last Sec- 
tion and write explicitly the Bogoliubov transformations 
in the amplitudes for the removal and addition transi- 
tions. For the removal amplitude this means. 
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(22) 



r 



An analogous expression may be written for the ad- ments of the operators 7 we have to adopt the following 
dition amplitude. To explicitly evaluate the matrix ele- approximations: 

I 



(GSA±2|7i,cri7«',o-2|GSA) ~ (GSA|7lo-i7n',<T2 |GSa) = ^rm"5cri<T27 



(GSyl±2|7 



IGS, 



(GSA|7n,<Ti7n 



(GSA±2|7i..,7l^.JGSA) ^ (GSA|7l.,7l^.JGSA) = 0, 



(23) 



and thus obtain for the removal and addition ampli- 
tudes 



J 



J-,^(2, + l)u;t-(rX(r), 



nlj 



J-j5:(2, + lX(r).„^,f(r). 



475-7-2 



nlj 



(24) 



The improvement with respect to the previous deriva- 
tion is achieved by actually making the approximation 
\GSa) ^ \GSa±2) on the matrix elements of the cre- 
ation/annihilation operators 7^,7. That is, the approxi- 
mation is done after having written and explicitly applied 
the Bogoliubov transformations in the expression of the 
transition amplitudes. This leads to expressions where the 
wave functions of the two nuclei appear as originally used 
in ref. !l] . 

The pair-transfer strengths are given by 



One notices also that 



r>Rem 



(A- 



pAdd 
^GS 



(A). 



(27) 



B. Canonical formalism 



pRcm 
^GS 



(A) 



dr^(2J + l)^;^,7^(r)^;4(r) 



nlj 



dr}_^(2j + l)<,^.(r)<+^(r) 

nlj 



1 (25) Similarly to what is done above in the quasiparticle 

formulation, a better approximation can be obtained also 
in the canonical formalism. Assuming that the canonical 
states do not change too much between the nucleus A 
and A ± 2 (no core polarization due to the addition or 
(26) removal of two nucleons) the improved formula are 



7 



and 



-^(GSA+2|aI4|GS^) (|^,(r,t)P + |v'.(r,;)l' 



i>0 



i>0 



(28) 



i>0 



Accordingly, the two-particle addition and removal prob- 
abilities now read 



J 



r 



(29) 
(30) 



It should be noted that the approximation made to ob- 
tain the above expressions slightly differs from the one 
used in the previous subsection and differences between 
the two sets of results may actually be expected within 
the improved formalism. 



C. Results and discussion 

Since the difference between removal and addition 
probabilities is just a shift of nuclei (see Eq. ([27)) '). we 
consider in what follows only removal probabilities. In 
Figs, m and [5] the results obtained with the improved 
expressions for the removal strength are presented for 
the mixed and the surface interaction, respectively. The 
canonical basis results (filled circles) are compared with 
the quasiparticle results (dashed line) and with the re- 
sults obtained with Eq. ^ (solid line) . In figure HI we 
see that the canonical and quasiparticle improved forms 
give similar results. In particular, these new forms of 
the pair-transfer probability lead to non-zero values for 
magic nuclei. It is also worth mentioning that the prob- 
ability is enhanced in the mid-shell with respect to the 
corresponding value obtained using Eq. ([B]) . 

The case of Figure [S] is different. We observe some ir- 
regularities in the trend of the results which have been 
obtained with Eq. ([25]) . These irregularities are related 
to some specific features of the HFB calculations in coor- 
dinate representation. It may happen that the fragmen- 



tation of the occupations among the discretized quasipar- 
ticle states (box boundary conditions with a box radius 
equal to 20 fm) is not the same in the nucleus A zL 2 
and in the nucleus A. When this is the case, some ir- 
regularities appear in the transfer strength calculated by 
making products of functions belonging to the different 
nuclei A and A ±2. This situation does not occur when 
one uses the canonical basis, that displays a smoother 
behavior when moving from one system to the neighbor 
one. There irregularities are much more pronounced in 
the case of surface-peaked interaction. These results in- 
dicate that the improved treatment of pair transfer dis- 
cussed here should be done using a discrete basis instead 
of the coordinate representation to avoid the irregulari- 
ties. 



We also observe that the differences between the results 
obtained with the canonical and the quasiparticle formu- 
lations are more pronounced in the case of a surface- 
peaked interaction (Fig. 5). This may be related to the 
fact that the artificial irregularities of the quasiparticle 
results are more important in the case of a surface pairing 
interaction (the fragmentation of the occupation numbers 
may vary more strongly from the nucleus A to the nucleus 
A ± 2 in this case) . The slightly different approximation 
adopted for the ground states of the two nuclei in the two 
derivations also accounts for these differences. 
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Mixed interaction 




Figure 4: Comparison of the removal probability obtained 
with the mixed pairing case using Eq. ([6]) (solid line) and the 
improved expressions given by Eqs. (|25p (dashed curve) and 
(filled circles). 



Surface interaction 
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Figure 5: Same as in Fig. |4] but for the case of a surface- 
peaked pairing interaction. 



V. EFFECT OF PARTICLE NUMBER 
CONSERVATION ON PAIR TRANSFER 
PROBABILITY 



The quasiparticle states which are generally used to 
estimate pair-transfer probabilities are not eigenstates 
of the particle number operator. In particular, a state 
IGSyi) contains not only components with the correct 
number of particles, but also components with particle 
numbers A±2, A±A, ... inducing spurious contributions 
to the probabilities. These contributions can be exactly 
removed using projection techniques. In the following, it 
is assumed that the ground state of a system with particle 
number A is written as 



\GSa)^\A)=P^I[[ 



i>0 



) |0), (31) 



where denotes the projector onto the number of par- 
ticles A. The projected state has the particularity that 
it has the same canonical basis as the original quasipar- 
ticle vacuum from which it is constructed. As shown 
in Appendix \^ the transition probability accounting for 
particle number projection can now be approximately es- 
timated by replacing occupation numbers of the quasi- 
particle state nf in Eqs. and by the new equiv- 
alent occupation numbers, denoted by nf, in the pro- 
jected state, i.e.: 



pAdd 

^GS 



pRem 



(32) 
(33) 



In Figs, inland [71 the effect of particle number conser- 
vation on the estimation of the transition probabilities 
is illustrated for mixed and pure surface pairing, respec- 
tively. To calculate the occupation number of the pro- 
jected state, standard gauge angle integration technique 
has been used with a Fomenko \1M discretization using 
199 points (see for instance Ref. jl4|). 

In these figures, we clearly see a very interesting and 
unexpected effect: at mid-shells, the projection actually 
tends to cancel out the effect of using the occupation 
numbers of the A and A ±2 nuclei in the improved for- 
mulation of last Section. The pairing strength gets again 
closer to the_probability obtained with the approxima- 
tion of Ref. is] and of Sec. III. However, close to magic- 
ity, projection has little effect and the results remain un- 
changed with respect to the improved treatment of the 
transfer probability. We mention that in Ref. [9], in 
order to have a non zero probability at shell closures, 
the transfer probability is calculated in these cases with 
the particle-particle random-phase approximation. Non- 
zero values are obtained here within a different approach 
based on a unified model suited to treat all the nuclei, 
both at shell closures and in mid-shell regions. However, 
the strength at shell closure is much lower here compared 
to Ref. 



VI. CONCLUSIONS 

Starting from an approximate derivation of the pair- 
transfer probability within the HFB approach, we have 
investigated its range of validity. Two different pairing 
interactions with a different surface/ volume mixing are 
used. To get a deeper and complementary physical in- 
sight we have formulated the same problem also in the 
canonical basis representation. In this case, the radial de- 
pendence of the wave functions is completely integrated 
out and disappears in the transfer probability. The same 
surface effects which are due to the integrated radial pro- 
files of the contributing wave functions in the quasiparti- 
cle representation still exist in the canonical basis case as 
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Mixed interaction 




Figure 6: Removal transfer probability obtained accounting 
for particle number conservation (open square) and compared 
to the improved (Sec. IV) (filled circles) and more approxi- 
mated (Sec. Ill) expressions (solid line) for the mixed inter- 
action. 
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Figure 7: Same as figure [6] for the case of pure surface pairing 
interaction. 

a consequence of the different fragmentation of the single- 
particle canonical occupancies. The transfer probability 
is actually connected to the diffuseness/ fragmentation 
of the single-particle state occupancies close to the Fermi 
energy. 

The possibility to improve the description of pair trans- 
fer probabilities in a mean-field model is then discussed. 
Two improvements are introduced: (i) the use of two 
different quasi-particle vacuum for the initial and final 
ground states; (ii) the possibility to exactly restore the 
proper number of particles in the entrance channel. It is 
shown that these refinements are important especially 
close to the magicity where a non-zero probability is 
found contrary to the simplest method. It turns out 
that the enhancement observed in the mid-shell by tak- 
ing two different ground states is partially compensated 
by the particle number restoration. We can thus finally 



conclude that the simple formula used in Ref. [9|, al- 
though approximated, can still provide good results at 
mid-shells. We remind however that the projection is 
performed here a posteriori. If variation after projection 
is made, additional correlations are expected to appear 
especially close to magic numbers [T7| . This open new 
perspectives and could be an interesting subject for fu- 
ture investigations. 
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Appendix A: Transfer probability for states with 
good particles number 

In this section, we consider two states denoted by \A) 
and 1^ + 2) obtained by projecting quasi-particle states 
onto good particle numbers, i.e. 

\A)^ ^ P^l[[l + 4alal)\0) (A2) 

i 

\A + 2)a+, ^ P^+^ll(l + 4+'alaD\0) (A3) 



where, we have introduced the coefficients xf — vf/uf 
and = vf^'^/uf^'^. The addition transfer prob- 

^ Note that, consistently with the approximation made in section 
IIII Bl it is impUcitly assumed that this many-body states shares 
the same canonical single-particle basis. 

^ Here, the convention \A) g is taken for the state: 



\A)b = P^Yl(l + xf 44)10) 



(Al) 



ability accounting for the particle number conservation 
requires to estimate the quantity a+2{A + 2\alat\A)A- 
An approximate form is derived below. 



1. Some properties of projected states 



The projection has the effect to select the component 
with good particle number in the quasiparticle vacuum 
leading to: 

where the B refers to the fact that the quasi-particle vacuum has 
been obtained with the constraint that the number of particle is 
fixed to B in average. 



\A)a ^ 
\A + 2)a+2 = 



N+1 



|0) 



(A4) 
(A5) 



where N and N + 1 denotes the number of pairs. The properties of projected state have been recently reviewed in 
Refs. [l^ [l^ and some of them will be recalled below. For instance, it has been shown that several recurrence relation 
exists to manipulate this state. Indeed, by developing the power in previous expressions and noting that (a|at)^ — 0, 
we deduce: 



N 



N-1 \ 



\A)a = 



}\0) 



= \A : i)a + VNxfajallA ~2:i)A 
from which we obtain the recurrence relation on the overlaps: 

{A : i\A : i)A = {A\A}a - N\xf\^{A ~2:i\A-2: i)a. 
By definition, all the states |. : «)b do not contain the pairs a\a\. Accordingly, we have: 



a\a\\A)A = \A:i), 



A+2{A + 2\4a\ = ^{N+l)xt+^*A+2{A:i\. 
With these expression, we finally see that we have the relationship: 



A+2 



{A + 2\a\4\A)A = ./W+^)4^^*A+2{A:i\A:t)A. 



(A6) 



(A7) 



(A8) 
(A9) 



(AlO) 
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2. Approximate form of the transition probability 

States defined by Eqs. (jA2l - IA3p are non-normalized states. To estimate transition density, one should use the 
following normalized states instead: 

leading to: 



A+2{A + 2\ala^A)A = ^{N + l)xf+^* ^ {A : i\A : i) a- 

^{A\A}a{A + 2\A + 2)a+2 



Assuming that 



A+2{A : i\A : i)A ^ A+2{A : i\A : i)A+2 ^ a{A : i\A : i)A (A13) 
and using the Eq. ()A7p . two approximate forms of the transition can be obtained: 

A+2(Ar2\a\a\\A)A ^ -j=^^^Z^a+2{A: ^\A■. ^)A+2 

^{A\A)a{A + 2\A + 2)a+2 

and 

A+2{Ar2W4\A)A ^ ^=^^^^=l^{{A\A)A~N\xt\^A~2■.^\A-2■.^)A). 

y'\A\A)A{A + 2\A + 2)^^2 

Combining these two expressions, it could be checked that: 

,t^tiJ\.|2 ^ , 1M^A+2.|2 (^:»I^:»)a+2 ^, ( , ^a,^ {A - 2 : i\A - 2 : i) 



\A+2{A + 2\ala]\A)A\' ^ + l)kr^1'7^-T^^r-^F^ x U - ^Nfl 



{A + 2\A + 2)a+2 V {A\A)a 

We recognize in this expression nothing but the occupation probabilities of the state i in the projected state respectively 
with A and A-l-2 particles [11] given by: 



and 



n. 



Altogether, we obtain that the transition amplitude between two states with good particle number can be approxi- 
mated by 



A+2 



{A + 2\a\a]\A)A ^ J nf+\l - nf) (A16) 



that is nothing but the same expression as the one obtained in the case of non-projected state except that the 
occupation numbers entering here are those associated with the projected states. 
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